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ABSTRACT 

To investigate the structure and composition of the dusty interstehar medium (ISM) of low surface 
brightness (LSB) disk galaxies, we have used multiwavelength photometry to construct spectral energy 
distributions for three low-mass, edge-on LSB galaxies (K-ot= 88 to 105 km s~^). We use Monte Carlo 
radiation transfer codes that include the effects of transiently heated small grains and polycyclic 
aromatic hydrocarbon molecules to model and interpret the data. We find that unlike the high 
surface brightness galaxies previously modeled, the dust disks appear to have scale heights equal to 
or exceeding their stellar scale heights. This result supports the findings of previous studies that low 
mass disk galaxies have dust scale heights comparable to their stellar scale heights and suggests that 
the cold ISM of low mass, LSB disk galaxies may be stable against fragmentation and gravitational 
collapse. This may help to explain the lack of observed dust lanes in edge-on LSB galaxies and 
their low current star formation rates. Dust masses are found in the range 1.16 to 2.38 x IO^Mq, 
corresponding to face-on (edge-on), F-band, optical depths 0.034 <; Tface ^ 0.106 (0.69 ^ Teq ^ 1.99). 

Subject headings: galaxies: ISM — galaxies: spiral — radiative transfer 
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1. INTRODUCTION 

Low surface brightness (LSB) disk galaxies, whose cen- 
tral, face-on surface brightnesses fall below fj,B,o ^ 23.0 
mag arcsec"^ in the B band, account for an important 
fraction of the lumin osity and galactic mass densities 



mass stars w ill slow metal production and indeed Mc 



of the local universe ( Driver 



1999 Minchin et al. 2004 



O'Neil fc Bothun|2000p . LHB dis k galaxies show a range 



of masses and mo rphologies (e.g., Sprayberry et al. 1997 



Auld et al.| 2006), but the most common are bulgeless, 
late-type LSB disks with Hubble types Sd-Sdm. These 
galaxies frequently show signs of bein g dark m atter dom 
inated at nearly all radii (e.g., de Blok & McGaugh 1997 
de Blok fc Bosma||2002]|Baiierjee erai.||201U^ 



LSB galaxies have been extensively studied in the Hi 
21-cm line. Often their total Hi masses are compara- 
ble to or larger than their high surface brightness (HSB) 
counterparts, while their H i surface densities are lower 
by a factor of ~2 (van der Hulst et al. 1993). In fact 



LSB galaxies are usually found to have most of the ir Hi 



disk at surface densities below the Kennicutt ( 1989 ) crit 
ical values for th e formation of massive stars! van der 
Hulst et al. 1993). This fact may lead to a lower than 



expected star formation efficiency causing LSB galaxies 
to evolve more slowly than HSB galaxies. A lack of high 
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Gaugh (1994) find LSB galaxies have metal abundances 
ot ~ 1/3 solar. Hii regions are o bserved in LSB galaxy 



1996 McGaugh 



disks as traced by Ha emission ( Rand] J 

et al. 1995 Matthews et al.|[l999 ) and blue dis k colours 
indicating young stars are present (McGaugh et al.|1995 



Matthews et al _ 

LSB galaxies do have ongoing star formation 



1999[ [Matthews fc lJson|2008| ) suggesting 

However, 



the conditions and structure of the interstellar medium 
(ISM) in LSB disk galaxies are still uncertain. 

While observations of Ha and H I have revealed infor- 
mation on the ionized and atomic gas, little is known 
about the dust and molecular hydrogen in the ISM. 
Based on observat ions of ^^CO(l-O ) emission in several 
LSB disk galaxies [Matthews et al.| ([2005) find that the 
CO emission, and hence the molecular hydrogen content, 
depends strongly on the rotational velocity. It appears 
that the rotational velocity, as well as the surface den- 
sity, play an important role in shaping the structure and 
co ntent of the ISM. This i s also suggested by the work 
of Dalcanton et al. (2004) who find that the gas and 
dust disk ot a galaxy is more stable to the effects of ra- 
dial instabilities at lower rotational velocity. This limits 
the ability of perturbations to cause the collapse of the 
ISM into a thinner layer and may inhibit global star for- 
mation by preventing the formation of large numbers of 
giant molecular clouds. 

Due to the difficulties of directly detecting molecu- 
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lar gas in low-mass LSB spirals (cf. [Matthews fc Gao 



2001 Matthews et al. 



tive method to inves' 



2005 Das et al. 



an alterna- 



igate the composition and struc- 
ture of their molecular ISM is to observe far-infrared 
(FIR) emission produced by dust. The dust and gas 
are expected to be well mixed and so the highest dust 
densities should trace the high density, cool phases of 
the ISM. In order to interpret the distribution of FIR 
emission we require a radiation transfer model to pro- 
duce synthetic images and spectral energy distributions 
based on galactic stellar and dust distributions. Pre- 
vious modeling efforts have focussed on fitting the op- 
tical aiid_jiear;infrared (NIR) extinc tion of HSB galax- 
ies jKylafelTB ahcall 1987; Xilouri s et al||1997[ ]l998 
1999| [iJiancETet al.||2000[ |Bianchi| |2007| ) and treatmg 
the dust absorption and reemission processes to model 
the full UV/optica l to sub-m ni spectral energy distribu- 
tion (S EP) (Popescu et al.l igOOO; Misiriotis et al...200ri 
Bianchi 2008; Bacs ct al. 20TU{ ]Popescu et ai.||20T'T| ~Tn 
the present work, we cxtend~tnis type of analysis tor the 
first time to a sample of LSB galaxies. 

When modeling the optical appearance of a disk galaxy 
an edge-on orientation is advantageous. This viewing ge- 
ometry allows both the vertical and radial ISM structure 
to be constrained simultaneously and small scale features 
of the galaxy are smoothed out due to overlapping sight 
lines. The effects of dust become more apparent and 
dark dust lanes can often be seen obscuring stellar emis- 
sion where the central regions of the galactic disk have 
become optically thick. This makes the quantification 
of dust masses and distributions much more straightfor- 
ward than the more degenerate case of a less inclined 
disk that is optically thin across most of its area. 

In this paper we further constrain the effects of dust, 
ISM geometry and mass on the UV/optical - FIR SEDs of 
edge-on LSB galaxies using Monte Carlo radiation trans- 
fer techniques. In section [2] we will introduce the sample 
of galaxies to which the analysis is applied followed by 
a description of the the data available at various wave- 
lengths in section [3] Section [4] includes a description of 
the model used in the analysis. The results for our sam- 
ple of galaxies can be found in section [5] and a discussion 
of the key results is presented in section |6] 

2. GALAXY SAMPLE 

We have selected three LSB galaxies that have been ob- 
served at a wide range of wavelengths to allow the best 
determination of the shape of the SED. The galaxies are 
nearby (<10 Mpc) and are well-resolved in optical imag- 
ing . We have selected the galaxies to span a range of 
central face-on surface brightness {^Bfl = 22.6—23.6 mag 
arcsec"^) as well as rotational velocity (Kot = 88 — 105 
km s~^). The three galaxies have similar total H i masses 
but observations indicate that their molecular hydro- 
gen co ntent decreases with de creasing rotational velocity 
(Frot) ( [Matthews et al.||2005[ ), suggesting that the ISM 
conditions may be ditterent across our sample. Table [T] 
contains the properties of the galaxies as gathered from 
the literature. High resolution optical imaging of the 
three galaxies can be seen in Figure [l] 

2.1. UGC 7321 



UGC 7321 is one of the best studied edge-on LSB 
galaxies. It shows a very thin disk (a/b ^ 10.3jr] with 
signs of clumping of the stellar and dust dis tributions in 
the inner galaxy ( Matthews fc Wood] 2001) and a more 
uniform structure m the outer disk. Tne disk shows 
strong {B — R) colour gradients in both horizontal and 
vertical directions that cannot be e xplained by the effects 
of dust (Matthews fc Wood[2001 1 and appear to be due 



to stellar population or metallicity gradients within the 
disk itself. Ha imaging shows ionized gas is present in a 
very thin layer that has very little structure and shows no 



signs of large scale star forming regions ( Matthews et al. 
1999|). The Hi gas content of UGC 7321 has a greater 



radial extent th an the optical disk of the galaxy ( Uson fc 
Matthews 2003 ) and sho ws indications tha t the Hi disk is 



also vertically extended (Matthews fc Wood 2003). Ob 
servations of the CO(l-0j line emission by Matthews 
Gao ( 2001[ ) confirm the presence of molecular gas in U G 
7321 and this measurement leads to a molecular hydro 
gen mass of Mh^ ^ 3.2 x IO'^Mq assuming a Galac 



tic CO-H2 conversion factor X = 3.0 x 10 
kms~i)-i( |Young fc Scoviiie|[l99l[ ). 



-20 



cm 



(K 



2.2. IC 2233 



IC 2233 is another example of a thin, bulge-less galaxy 
viewed close to edge on. It shows no signs of a dust 
lane in the central regions of its disk with only a few 
dark clouds visible. Rand ( [T996j ) detected Ha emission 
associated with Hii regions across the disk of IC 2233, in- 
cluding a large complex in the outer region of the galaxy, 
indicating that there is ongoing star formation. In con- 
trast to UGC 7321 the H 11 regions in IC 2233 are often 
located out of the mid-p lane at up to ~ 0.5Zs i n the 
outer parts of the galaxy ([Matthews fc Uson[[2008| . Ob- 
servations by Matthews fc Uson" pnnS| find that Hi in 
IC 2233 is extended vertically and horizontally beyond 
the optical disk, with evidence for flar i ng in the outer 
gas disk. Attempts by [Matthews et al. (2005) to detect 
CO emission from IC 2233 were unsuccessful leading to 
a 3-cr upper limit on the mass of molecular hydrogen of 



1.4 X lO^Mp 



2.3. NGC 4244 



NGC 4244 is the closest galaxy in our s ample lying 
at an estimated distance of only ~ 4.4Mpc (Seth et al. 
2005a). NGC 4244 shows many dark clouds in optical 
imaging of its c entral regions. Despite this, based on 



NIR photometry Kodaira fc Yamashita ( 1996 ) conclude 
that the disk of NGC 4244 sutters from little internal 
ex tinction. Continuum subtracted Ha imaging obtained 
by Hoopes et al. ( 1999 ) shows H 11 regions extending up 
to ~ 990pc from the disk midplanq^ The two largest 
H II regio ns are found in the outer galactic disk. The 
work of Oiling (1996) investigating the Hi distribution 
in NGC 4244 shows that the gas distribution appears 
to be slightl y extended radially comp ared to the stellar 
disk. Both (Matthews fc Gao] ( |2001[ ) and [Sage] |1993 1 
report the detection ot CO emission from NGC 4244 and 
the estimated H2 mass is 1.4 x 1O''M0. 

^ a/b is the disk axial ratio measured at the 25.0 mag arcsec"'^ 
i?-band isophote 

■^ This value has been re-scaled to the distance adopted in this 
paper (4.4Mpc) 
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IC 2233 have previously been observed using the Wide 
Field Camera (WFC) of the Advanc ed Camera for Sur - 



Figure 1. High resolution imaging of the inner regions (approx- 
imately ^ 3.4 X l.Skpc) of the three galaxies studied. WFPC2 
F702W and FSMW composite image of UGC 7321 (top). ACS 
WFC F606W filter images of IC 2233 {middle) and NGC 4244 
(bottom). The white bar in the top right corner of each image has 
a length of 500pc, assuming the distances in Table IT] UGC 7321 
and NGC 4244 show the presence of dark clouds in their nuclear 
regions that are seen in lower numbers in IC 2233. This may be 
indicative of a larger mass of molecular gas in these galaxies which 
could have an impact on the FIR emission if dust is also associated 
with these clouds. 

3. DATA 

All three galaxies have been observed in the optical 
by the Sloan Digital Sky Survey (SDSS) in u, g, r, i and 
2 band s and the data are ava ilable under data release 
severr (Abazajian et al 



loadea 



2009 1 . The images were down- 
rom the SDSS archive and fluxes in each band 



were extracted using GAIA's aperture photometry tool. 
It was necessary to re-analyse the SDSS images to extract 
reliable photometry as each galaxy appeared to have 
been split into multiple components during the SDSS 
photometric reduction. 

We have also made use of B band imaging of UGC 
7321 and IC 2233 obtained at the 3.5m WIYN telescope 
at Kitt P eak, AZ. The se observations are described in 
Matthew s_et al.| ( |1999[ ) (UGC 7321) and [Matthews fc 
Uson, (,20d8f ■■■'"^—^ — ■ ' ^' " - 



ibratC' 



(1C2233). The images have been hux cal- 
and are of a higher resolution than the SDSS 
imaging available. In the case of NGC 4244 we do not 
have access to higher resolution data and so the SDSS r 
band imaging will be used for the optical comparisons. 

High resolution Hubble Space Telescope {HST ) imag- 
ing is also available for our three galaxies. NGC 4244 and 

4 http://www.SDSS.org/dr7/ 



veys (ACS) with the F606W filter (|Seth et al.||2005b[ ) 
These data were accessed from the Multimission Archive 
at STScI (MASTfl Imaging of UGC 7321 is available 
from the Wide Field Planetary Camera 2 (WFPC2) and 
we use a composi te F702W and FSl^W im age previ- 
ously presented in [Matthews fc Wood] ( |2001[ ). HST im- 
ages of the inner regions of each galaxy can be seen in 
Figure [T 

The MIR and FIR observations of our three galaxies 
were obtained with the Spitzer Space Telescope's IRAC 
and MIPS imagers. Data were obtained for this project 
under Spitzer PID 20432 (PI: L. Matthews) along with 
observations taken from PIDs 3 (PI: G. Fazio) and 40204 
(PI:R. Kennicutt). These data have recently been used 
as part of t he Spitzer Local Volume Legacy (LVL) survey 
( Dale et al.^^2009) . This provides imaging and photome- 
try tor a large sample of galaxies within ~llMpc. Data 
are available for all four IRAC bands (3.6,4.5,5.8 and 
S.O^m) as well as three MIPS bands (24, 70 and IGO^im) 
with photometry provided for each band after foreground 
stars and background galaxies have been removed and 
aperture corrections applied. 

Photometry is also available based on near-infrared 
(NIR) 2MAS9j^ imaging using apertures matched to 
those used in the extraction of the IRAC and MIPS 
fluxes. Spitzer LVL survey imaging can be accessed via 
an online archivqj and NIR, IRAC and MIPS fluxes are 



available in Table 2 of Dale et al. 



Lee et al. 
far and near 



( |2009[ ). 

(|2011[ ) provide a catalogue of matched 
iv ot)servations for galaxies contained in 



the Spitzer LVL survey that have been observed with 
GALEX. The GALEX bands are centered at 1528 A (far- 
UV) and 227lA (near-UV) ( JMorrissey et al.||2005[ ) and 
fluxes have been extracted using the apertures obtained 
from the FIR Spitzer data for consistency. The UV data 
provide a constraint on the emission from the young stel- 
lar population which, due to the higher opacity at UV 
wavelengths, will play an important role in heating the 
diffuse dust. 

4. RADIATIVE TRANSFER MODEL 

In order to produce synthetic images and SEDs for 
comparison to the data we adopt a Monte Carlo radiation 
transfer modelling scheme. 

4.1. Emissivity and Dust Distributions 

The model allows energy packets emitted within a 
smooth stellar emissivity distribution to be tracked as 
they propagate through a dusty medium. For both the 
emissivity and dust we adopt a smooth "double expo- 
nential" distribution: 



p{zu,z) ex exp(— |z|/Z) exp(— nj/i?). 



(1) 



where p{'cu, z) is the density at a point in the galaxy 
specified by w, the cylindrical radius, and \z\ the height 
above or below the galactic plane. The scale lengths R 
and Z can be independently varied for both components 

^ http://archive.stsei.edu/index.html 

^ http://www.ipac.caltech.edu/2mass/ 

''' http://irsa.ipac.calteeh.edu/data/Spitzer/LVL/ 
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Table 1 

Galaxy Properties 



Vrc 



P^bM" 



M. 



UGC7321 


Sd 


1.1 


10 


IC2233 


Sd 


1.8 


10 


NGC4244 


Sd 


3.6 


4.4 



M 



H2 



Name Type B band Distance 

Luminosity 

(IO^Lq) (Mpc) (kms-i) (mag/"^) (IO'^Mq) (IO^Mq) 

105 



95 



23.6 


1.1 


3.2 


22.6 


1.1 


<0.1 


23.6 


1.3 


1.4 



^value deprojected to face-on 

as necessary. The stellar emissivity and dust density are 
discretised across a three dimensional spherical grid of 
cells to allow our Monte Carlo radiation transfer routines 
to be used. To reduce the number of free parameters in 
the model no attempt has been made to take account 
of clumpy substructure within either the emissivity or 
dust distributions, which may be important in reproduc- 
ing the optical properties of galactic disks (iMatthews & 
Wood|2001[|Bianchi et al.|2000[[PiCTini et al. 2004). 'i'his 
also allows the use ot a lower resolution grid of cells, re- 
ducing the memory requirements of each radiation trans- 
fer simulation. 

4.1.1. Emission Sources 

To begin the radiative transfer process it is necessary to 
sample an initial wavelength of emission for each photon 
packet. Wavelengths are chosen to reproduce the intrin- 
sic unattenuated stellar SEDs produced by the popula- 
tion of stars present in each galaxy. For this purpose a 
set of synthetic stellar SEDs was constructed using the 
GALE^ codes ( [Kotulla et ari [2009). These are shown 
in Figure [2j Each model initially has a constant SFR 
for several Gyrs which is followed by an exponentially 
decreasing burst of star formation with an e-folding time 
of 1 Gyr. Table [2] lists the values used within the model 
for each galaxy. In all cases a Salpeter initial mass func- 
tion with slope a = —2.35 is used and the metallicity of 
the gas is allowed to change with time to remain chemi- 
cally consistent. We also include the effects of continuum 
and line emission from gas within the model SEDs. 

In searching for appropriate stellar populations to 
model the observed stellar SED we have found that a 
wide range of ages can provide adequat e fits. This has 
previously been noted b y other authors ( Zackrisson et al. 



2005 Vorobyov et al. 



2009) who have suggested that 



based on additional information, such as Ha equivalent 
widths and oxygen abundance gradients, that old stellar 
populations are more probable in samples of blue LSB 
galaxies. The ages adopted here range from 8.6 — 9.4 



Gyrs which are in agreement with Jimenez et al. ( 1998 1 
who suggest ages larger than 7 Gyrs tor Ltiti galaxies. We 
note, however, that our conclusions on the dust masses 
and distributions are most sensitive to the shape of the 
intrinsic stellar SED and will not be significantly altered 
by the exact age, metallicity or star formation history of 
the underlying stellar population which produces it. 

The SFRs for UGC 7321 and IC 2233 have previ- 
ously been estimated based on the observed Ha emis- 
sion. SFRs of ~ 0.02 and 0.05 Mp,yr~ ^ have been found 



for UGC 7321 (Matthews & Wood 20031 and IC2233 
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Figure 2. The intrinsic stellar emission templates used for our 
three galaxies. The fluxes have been scaled for comparison. The 
important property of each template is its shape, as this is used 
within our models to correctly sample the stellar emission. Each 
galaxy required a different intrinsic stellar template to match the 
observation in the UV . optical and NIR. T he templates were pro- 
duced using GALEV jKotulla et al.|[2009[ l and model parameters 
can be found in Table 



* http://www.galev.org/ 



( Matthews fc Uson||2008 ) respectively, which are in rea- 
sonable agreement with the values in Table [2J We have 
estimated the SFR of NGC 4244 based on the Ha lumi- 
nosity, Lffa — 10'*''erg s~^ ( Kennicutt et al. 2008), and 
the SFR - Ha relation taken from Kennicutt| (il998|) . To- 
gether these yield a SFR - 0.079 ±0.023 M^yr-^ which 
is roughly consistent with the value adopted for the tem- 
plate stellar SED. 

A single stellar SED is assumed for each galaxy mod- 
eled. No account of stellar gradients, which may be i m- 
portant in LSB galaxies (Matthews & Wood 2001 1 is 
made in order to reduce the complexity ot the model. 
We do not believe that the stellar gradients will have 
a significant impact on the SED and derived dust disk 
properties. 

Also included is a secondary emission component from 
obscured star formation regions that may be necessary 
to account for the total FIR emission. As such emission 
takes place from compact regions that are below the res- 
olution of the spherical grid currently employed in the 
radiation transfer, a separate model has been created. 
The regions are approximated by a single central star 
surrounded by a shell of dust and gas. The shell has a V 
band optical depth of ry ~ 88 and extends from a dis- 
tance d — lO^^'pc out to 2.3pc. The size of the shell was 
chosen to best reproduce the shape of the o bserved FIR 
emission of star forming regions observed by |Chini et al.| 
(1986). The shell is exposed to internal heating from a 
stellar sou rce, with emiss ivity sampled from a model at- 
mosphere (Kurucz 1993) of temperature T = 36, OOOK. 
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Table 2 

Intrinsic Stellar Template Model Parameters 



Name 


Initial SFR 


Burst Mass'' 


Burst Time'" 


Age'^ 


SFR'' 




Mgyr-i 


Fraction 


109 yrs 


lO'' yrs 


Meyr-i 


UGC7321 


0.02 


0.1 


7 


8.6 


0.02 


IC2233 


0.005 


0.1 


8 


9.1 


0.04 


NGC4244 


0.03 


0.2 


7.5 


9.4 


0.04 



''Fraction of remaining gas converted to stars during burst 
''Age of galaxy at the time of the star formation burst 
'^Age of the galaxy adopted to fit the broadband fluxes 
''SFR at adopted age 



4.1.2. Dust Model 
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Figure 3. The template for the emission from compact sources 
between 10 and 2000/im. The solid line shows the model output 
SED as it varies with wavelength. The model comprises of a cen- 
tral star {T=36,000K) heating a constant density shell of gas and 
dust surrounding it. The shell extends from 10~* — 2.3pc and 
has a V band optical depth of Ty ^ 88. The total luminosity in- 
cluded in each cloud model is 2.5 X IO^Lq and each has a dust 
mass of 5 X IO^Mq. Data points represent t he median observed 
properties of 56 compact H II regions found by |Chini et al.| | |1986[ l 
at 12, 25, 60, 100 and 1300/.tm. Error bars show the range ot values 
observed at each wavelength. 



The shell is also illuminated from the exterior by a radia- 
tion field approximating the mean intensity in the plane 
of the galaxy model. Due to the size of the clouds, how- 
ever, the luminosity incident on the exterior is typically 
much smaller than the luminosity of the central source 
and hence diffuse illumination has little effect on the 
overall SED. The output SED can be seen in Figure [3] 
along with observed d ata points of 56 Galactic compact 
H 11 regions taken from Chini et aT] (|T986). A single SED 
is used to approximate all compact sources in the galaxy. 
The spatial distribution of compact sources is an expo- 
nential disk (Eq. II]) with Z — 0.1 kpc and R — 1.8 kpc 
in all models. This secondary emission component, from 
obscured star forming regions, has previously been found 
necessary to rep roduce the FIR emission in HSB galaxies 
Popescu et al. 12000; M isiriotis et al.|200T [Popescu et al. 

unt tor such en 



uffi r 



other methods to account tor sucn emission in- 
elude utilizing sub gri d resolution models for individual 
star forming regions ( [Jonsson et al. 2010) or the inclu- 
sion of adaptive grids to allow radiation transfer on the 
small size scales needed. Current adaptive grids models, 
however, lack t he resolution t o resolve individual star 
forming regions (Bianchi||2008). 



The dust content of the model galaxy is treated in 
a two component manner as described by [Wood et al.| 
( 2008 ) in order to take account of both large grains and 
small transiently heated grains. Large grains with sizes 
> 200A are assumed to have attained radiative equi- 
librium temperatures. In the case of these grains the 



Bjorkman & Wood (2001) method is used to update the 
temperature ot the cell and correctly sample the emissiv- 
ity for the reemitted photon packet when an interaction 
takes place. The advantage of this method lies in the way 
the frequency of the reemitted photon is sampled to cor- 
rect the frequency distribution of the photons previously 
emitted from the cell to the new cell temperature. The 
SED of the cell then relaxes to the equilibrium value, as 
the simulation progresses, without the need for iteration. 
The composition and size distribution of the large grains 
is taken from the work of Kim et al. (1994) who used a 
mixture of silicate and graphite grains to ht the observed 
extinction in the Milky Way. This treatment of large 
dust grains, without the extensions added to treat tran- 
sient heating of very small grains (VSG) and polycyclic 
aromatic hyd rocarbon (PAH) m olecules, has previously 
been used by Savoy et al. ( 2009 1 to investigate the dust 
conten t of early-type galax ies. 



The Bjorkman & Wood (2001) method cannot, how- 
ever, be applied to transiently heated VSGs and PAH 
molecules that are not at equilibrium temperatures. 
Emission from these grains is important in reproduc- 
ing the MIR emission. For photon packet interactions 
with grains that are sma l ler th an 200A the methods in- 
troduced in Wood et al. (2008) are used. The emission 
from PAH molecules and VS Gs is calculated from pre- 
computed emissivity files from Draine & Li (2007). The 
emissivity is assumed to vary only with the mean in- 
tensity of the radiation field and not with the spectral 
shape. The cutoff in the grain size at 200A is somewhat 
arbitrary but it is thought that grains a bove this will 
norm ally reach equilibrium temperatures (Draine & Li 
2007 1 . In order to correctly sample the emissivity ot the 
PAH molecules and VSGs a value for the mean intensity 
in the cell where the interaction took place is needed. 
This requires an iterative process where initially a value 
of J = 1 is use d, J being the value of the mean intensity 



relative to the Mathis et al. 



( 1983 ) value of the local in- 



terstellar radiation field (ISRF) of Jisrf = 2.17 x 10 
erg cm"'^ s^^. Consecutive iterations then use the mean 
intensity computed for each cell on the previous itera- 
tion to calculate the emissivity of the PAH molecules and 
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VSGs. Due to the small contribution of photon packets 
reprocessed by PAH molecules and VSGs to the total 
mean intensi ty it is found that typically three iterations 
are required ( Wood et al.||2008 ). 

As the ISM conditions of Ltia galaxies are likely differ- 
ent from those found in the Milky Way, it is possible that 
our treatment of the PAH/ VSG dust components may be 
inappropriate. It is known that the strength of the PAH 
emission relative to that from VSGs, estimated using 
the /iy(8/im)//i,(24/im) colour, decrease s with decreas- 
ing m etallicity (Engclbracht et al. 2005; Madden et al. 
2006 1 . This is likely due to either the delayed torm ation 
ot pAH molecul es in low metallicity environments (Gal- 
liano et al.|2008 ) or de struction/processing caused by the 
harde r radiation field ( Gordon et al.|2008 Madden et al. 
2006 1 . It is therefore possible that our treatment, which 
assumes Milky Way abundances, may overestimate PAH 
emission relative to that of the VSGs. The fraction of the 
total dust mass contributed by PAH/VSGs is also fixed 
to a value calculated for the Milky Way. This could affect 
the relative strength of the MIR and FIR emission. 

To model the appearance and SED each galaxy we can 
alter the following free parameters of the model: stellar 
scale length, stellar scale height, dust scale length, dust 
scale height, stellar luminosity, total dust mass, inclina- 
tion and the luminosity of the obscured star formation 
component. 

The overall fitting process was done manually to find 
the set of parameters that could provide adequate fits to 
both the optical and the UV-FIR SED. 

5. RESULTS 

Initial values for the fitting process were taken from 
sources in the literature. Matthews & Wood (2001) have 
previously investigated IJGC 7321 using a three dimen- 
sional Monte Carlo scattered light code and their pa- 
rameters for a smooth stellar and dus t distribution were 
adopte d as initial values in t his case. Matthews fc Usoii 
(2008) and Fry et al. (1999") have previously estimated 
the scale lengths ot 1C2233 and NGC 4244 respectively, 
by fitting an exponential function to the observed sur- 
face brightness profiles. While the functions used take 
no account of the effects of dust on the galaxy profiles 
they should provide reasonable initial estimates for the 
stellar distributions assum ing that dust effects are not 
severe. From the work of Xilouris et al. ( 1999 1 on high 
surface brightness, edge-on, galaxies we initially adopted 
the relation that the dust scale heights would be approx- 
imately half the stellar scale heights and the dust scale 
lengths around 1.4 times larger than the stellar values. 

One of the main difficulties in modeling the dust distri- 
butions of LSB galaxies is that their disks appear to be 
optically thin, even when viewed edge on. When mod- 
eling the appearance of HSB galaxies, the effects of the 
dust distribution can be quantified based on the presence 
of a dust lane along the galactic mid-plane. As none of 
our galaxy sample shows any sign of a dust lane the prob- 
lem becomes more degenerate. Initially we increased the 
dust mass in the adopted geometry until the effects of 
the dust became obvious and the optical profiles of the 
model galaxy no longer matched the shape of the ob- 
served data. It was found that the dust mass required 
to reprocess the stellar light and reproduce the observed 
FIR emission could not be located in an average HSB 



type dust disk without revealing its presence by fiatten- 
ing the optical profile in the central region of the galaxy. 

We find that the optical appearance and total FIR 
emission of the LSB galaxies can be reproduced by re- 
moving the constraint that the dust be more concen- 
trated in the galactic midplane than the stars. If the dust 
scale height is taken to be equal to the stellar scale height 
then it is also possible to fit the optical appearance of the 
galaxy adequately. In this case, however, a significantly 
larger dust mass can be accommodated within the dust 
disk without the appearance of a dust lane. In order to 
reproduce the observed FIR flux distribution it was also 
necessary to increase the dust scale length. When the 
dust is vertically extended into a disk with a larger scale 
height that is well mixed with the stellar population, the 
exact radial distribution of the dust becomes poorly con- 
strained from optical data alone. We have found the dust 
scale length is between 1.8 — 2.6 times the stellar scale 
length for our LSB disk galaxies. 

The best fit parameters for our models of each galaxy 
can be found in Table [3] Rg and Rd are the radial scale 
lengths of the stars and dust while Z^ and Z^ are the 
vertical scale heights. Lg is the intrinsic bolometric lumi- 
nosity of the stellar population. Md represents the total 
dust mass and i the inclination angle. Ldoud gives the 
luminosity emitted fr om th e template of compact dust 
emission (see section |4.1.1 ) . Tjace and Teq give the V 
band face-on and equatorial optical depths respectiveljr] 

Figures l4J [s] and p] show a comparison of the B (or r)^ 
3.6/im , 8/im7YO/im and 160/^m data (left panel) and syn- 
thetic images (right panel). In general the B (or r) and 
3.6/im model images match the large scale morphology 
of the data reasonably well. Differences are, however, 
apparent in the longer wavelength emission at 8, 70 and 
160/im that is dominated by emission from PAHs and 
dust grains. IC 2233 and NGC 4244 show significant 
structure in their 8 and 70/im images, suggesting that 
perhaps a more complex non-axisymmetric dust distri- 
bution is needed. In both galaxies several small point 
like sources in the outer disk are prominent sources at 
70/im. These are likely star forming regions as they ap- 
pear brighter at shorter wavelengths and show associated 
emission in GALEX iai-\]V and Ha imaging. The 160/im 
emission from all three of our model galaxies shows a sig- 
nificant deviation from the data. The models predict a 
centrally concentrated 16G/im image while the data sug- 
gest a more diffuse distribution. Figures [7j [8] and [9] show 
B (or r) band intensity profiles parallel to the minor axes 



of the galaxies at various points while figures 
[12] show the output model SEDs for the three 
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Overall the models are able to reproduce the global 
properties of the data. The surface brightness slices par- 
allel to the minor axis and along the major axis show a 
similarity to the data and deviations are likely caused by 
the comparison of smooth axisymmetric models to obser- 
vations of a galaxy that shows a clear clumpy structure 
(see Figure [l| . The main discrepancies found between 



model grid perpendicular to the plane of the disk. Teq is measured 
from the disk center to the edge of the model grid parallel to the 
disk plane. 
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Figure 4. Observed (loft) and synthetic (right) images of UGC 7321 at, from top to bottom: B band, 3.6/^m, 8^m, 70/jm and 160/im. 





Figure 5. Observed (left) and synthetic (right) images of IC 2233 at, from top to bottom: B band, 3.6/jm, 8/jm, 70/^m and 160/^m. 





Figure 6. Observed (left) and synthetic (right) images of NGC 4244 at, from top to bottom: r band, 3.6/im band, 8/xm, 70^m and 160^m. 



the surface brightness of the model and data generally 
occur in the central regions of the galaxies. As can be 
seen in Figure [T] these are also the regions that show the 
greatest number of dark clouds and bright star clusters. 
Our smooth axisymmetric model is not able to reproduce 
such structures and so only a relatively poor fit can be 
achieved in the central regions. The short dashed line 
in Figure [lO] shows the result of fitting the optical imag- 



ing of UGC 7321 using t he average HSB gala, xy stellar- 
to-dust scaling relation of Xilouris et al. ( 1999 ) . In order 
to produce the required 70^m emission a large luminosity 
is assig ned to the obscured star formation template (see 
section 4.1.1 ). However, this scenario underestimates the 
160/ini emission by a factor ~ 2. 

In the SEDs the shorter wavelength data from the UV 
through to NIR originate in the underlying stellar popu- 
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Lcloud 


Tface 


Teg 


UGC7321 
IC2233 

NGC4244 


1.8 
1.7 
1.9 


0.12 
0.18 
0.2 


4.30 
4.50 
3.4 


0.12 
0.18 
0.2 


14.5 
19.3 
33.5 


2.25 
1.16 
2.38 


88.5 
88.5 
84.5 


0.67 
0.01 
1.0 


0.072 
0.034 
0.106 


1.99 
0.69 
1.57 
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Figure 7. The top two plots show UGC 7321 B band surface 
brightness slices parallel to the minor axis at the galaxy centre 
(top), and 3 kpc radial distance along the major axis (middle). 
The lower panel shows the surface brightness variation along the 
major axis. All slices are average values taken across a ^ 60pc 
region. Solid lines show the predicted intensity based on our best 
fitting Monte Carlo radiation transfer models and crosses represent 
the data taken from the B band image. 

lations and can be highly attenuated by the diffuse dust 
in the ISM. Most of the photons absorbed by the dust, 
in our models, originated as short wavelength UV emis- 
sion. This can be seen in Figure [TO] as the difference 
between the intrinsic stellar emission (long dashed line) 
and the observed emission (solid line) once the photons 
have propagated through the dusty ISM. The amount of 



120 

100 

80 

60 

40 

20 


50 



40 ■ 



30 ■ 




Parallel to minor axis 



Parallel to minor axis 
r = 3kpc 





-4-2 2 4 
Offset (kpc) 



Figure 8. As in Fig. [t] but for IC 2233. 



absorption tends to decrease with increasing wavelength 
as the dust opacity decreases and the photons are able to 
travel through the ISM with a lower probability of inter- 
action. Our models can reproduce the observed optical 
and NIR fluxes well but show slight discrepancies in the 
FUV emission. The models are also unable to reproduce 
the z band flux (0.893/.tm). We believe this may be due to 
the lower S/N in the z band combined with the intrinsic 
low surface brightness of the galaxies, leading to galaxy 
flux being lost during removal of the sky background. 
The 5.8 and 8.0/im IRAC bands trace emission from 
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Figure 9. As in Fig. [71 but for NGC 4244 and utilizing SDSS r 
band rather than B band data. 

PAH molecules while the MIPS 24/ini emission is pro- 
duced primarily by the warm VSGs. In all three cases 
the 8.0/im emission is over-predicted by our models. The 
24/im flux is over-predicted by the model of UGC 7321, 
under-predicted for IC 2233. This behavior in the MIR is 
likely a result of our treatment of the PAH/VSG emission 
which is based on Milky Way abundances and illuminat- 
ing radiation field shape (see 4.1.2). Both the relative 
abundances of PAH molecules and VSGs in LSBs and 
the shape of the ISRF in LSBs are likely different caus- 
ing the discrepancies observed in the MIR. 

The larger, cooler dust grains are responsible for the 
peak of the FIR emission bracketed by the MIPS ob- 
servations at 70 and 160/xm. In all three cases we have 
been able to reproduce the FIR 70 and 160/xm emission, 
within the quoted photometry errors, using our models. 

6. DISCUSSION 

The model outputs suggest that it is possible to provide 
a good global fit to the optical appearance of LSB galax- 
ies, while at the same time reproducing the observed FIR 
emission. At some locations, however, such as the cen- 



tral region of UGC 7321 (Figure [7J top) there is a clear 
discrepancy between the data ana model. In most cases 
this occurs due to structure in the galaxy that is not in- 
cluded in our models, such as star clusters, Hii regions, 
and optically thick dust clumps. 

6.1. Implications of the Models 

In our analysis we have found that in order to accom- 
modate the mass of dust required to reprocess the FIR 
emission seen in our sample of three LSB galaxies it is 
necessary for the dust to be located in a diffuse disk 
with a scale height at least equal to that of the stellar 
population. This is in contrast to the dust distributions 
obtained by previous studies of HSB galaxies (see Section 
[5]) , which find the dust to be located in a diffuse disk with 
approximately half the scale height of the stellar disk. 
The main observational signature of the smaller vertical 
extent of the dust in HSB galaxies is the central dust 
lane observed when the inclination approaches edge-on. 
Dalcanton et al. (2004) have previously noted that in a 
sample of edge on galaxies that spans a range in Vi-ot there 
appears to be a remarkably sharp divide between galaxies 
that contain dust lanes and those that do not at around 
Kot — 120kms~^, with no abrupt change in the stellar 
or H I surface densities. They attribute this sharp divide 
to the onset of radial instabilities in galactic disks caus- 
ing the vertical collapse of the ISM when perturbed by 
spiral structure. The collapse leads to a higher gas den- 
sity in the central galaxy having the effect of increasing 
the star formation rate and surface brightness. The in- 
creased gas density also causes a corresponding increase 
in the dust density in the galactic mid-plane which, when 
viewed edge-on, will produce the pronounced dust lanes 
observed. From Table [l] it can be seen that all three 
of our LSB g alaxies have Vj-nt < 120kms~^ and so from 



the results of Dalcanton et al. ( 2004 ) would be expected 



to be stable across most of their disks. This result sug- 
gests that as the disks are stable the gas and dust disk 
should have a scale height comparable to or larger than 
the stellar value, which is what we have found from our 
modeling of the optical, NIR and FIR data. 

In the case of the LSB galaxies studied here an al- 
ternative explanation for the comparable scale heights 
of the dust and stellar disks is that the scale height of 
the stellar disk is comparatively smaller than that found 
in HSB galaxies, contributing to the thin appearance of 
the galaxies. In this case the dust and gas disk scale 
height may be comparable to more massive galaxies but 
a drop in stellar scale height has concen t rated the stars 
within the dusty ISM. Bizyaev & Kajsin (2004) have ob- 
served that stellar disks appear to become thinner as the 
R band central surface brightness dims, which provides 
evidence that we may indeed be seeing a reduced stellar 
scale height in the LSB galaxies studied here. This ef- 
fect may be due to the inefficiency of processes that con- 
tribute to the vertical heating of the stellar component 
such as scattering by giant molecular clouds, the effects 
of bar formation and a lack of gravitational interactions 
with neighboring galaxies. 

Alternatively we may be observing a combination of 
both the processes mentioned above. An increase in the 
dust disk scale height caused by the stability of the ISM, 
coupled with a relative flattening of the stellar disk. Both 
effects combined would then result in comparable scale 
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Figure 10. The observed and modeled SED of UGC 7321. Square points indicate an observed flux from one of our data sources either 
GALEX, SDSS, 2MASS or SpitzerlHAC/MlPS along with their associated errors. Error bars associated with SDSS imaging do not include 
the contribution of the sky subtraction errors. The solid line shows the model output SED as calculated from the best-fit to the WIYN B 
band image. Open circles show the predicted fluxes for Spitzer IRAC/MIPS bands. The long dashed line indicates the input, unattenuated, 
stellar template plus compact dust emission templa te. The dotted line shows the SED from the best-flt model when an average stellar-dust 
distribution for HSB galaxies l|Xilouris et al.|1999[l is used to flt the optical imaging. 



10 



-10 



E 

0) 



10 



-11 




0.1 



1 10 100 

wavelength {\im) 



1000 



Figure 11. The above plot shows the model SED of IC 2233 (solid line) along with the available data (squares). Open circles indicate 
the predicted model fluxes for the filters used in the observations. The dashed line shows the intrinsic SED of the emissivity sources (both 
stellar and dust emission from unresolved sources). 
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Figure 12. Same as Figure [TT] but for NGC 4244. 
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heights of the dust and stellar distributions. 

Previous works have suggested that the dust scale 
length in galaxies may range from being only s lightly 
larger than the stellar scale length (Muoz-Matcos et al. 
2009 1 to being mod erately extended ( piTlouris e t ai. 19991 
I'empei et al.|2010| or even u p to an order o± magnitude 
larger (Holwcrda ct al.||2005 ). The average value of dust 
scale length we hnd is K^ = 2.3±0.24 Rg which is slightly 
larger than previous values based on FIR and optical ex- 
tinction modeling b ut significan t ly sm aller than the value 



7321 and NGC 4244 faU between those found for the 



al. (2005) from background 



estimated by Holwerda et 
galaxy counts. 

The low star formation rates observed in LSB galaxies 
are likely caused by the vertically and radially extended 
gas distributions leading to low surface and volume densi- 
ties that inhibit star formation via the formation of cold, 
dense clouds 



van der Hulst et al. 
peak Hi surface density, in a sample 'o 
ies, of4.6±1.4M,-— -2 
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1993) find a mean 

SB disk galax- 

This value is around a factor ^ 2 



lower than that found by Cayatte et al. 
Although the low density ot' 
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galaxies. 

leads to gravitational stability and prevents global star 
formation driven by spiral density waves, star formation 
still occurs in localized high density regions. Ha imaging 
suggests that this does not occur in the same pattern in 
all LSB disk galaxies. U GC 7321 shows Ha emi ssion con- 
centrated in a thin layer ( Matthews et al.|1999 l) while IC 
2233 and NGC 4244 show more vertically extended emis- 
sion and large star forming complexes in th e outer disk 
([Matthews fc Uson|2008] [Hoopes et al.|1999| ). These dif- 
lerences in the morphology ot the current star formation 
suggest that there may be differences in the clumpy ISM 
structure. This may indicate a more gradual transition 
in the ISM structure across galaxies with decreasing Vrot 
rather than simply an abrupt change at a given galactic 
mass. 

It can also be seen from Tables [3] and [T] that in order 
to model the FIR emission of UGC 7321 and NGC 4244, 
which are known to contain molecular hydrogen in the 
nuclear regions of their disks, it has been necessary to 
include a larger fraction of emission from our template 
of embedded compact star forming regions. IC 2233 on 
the other hand has not been detected in CO emission 
and, indeed, is shown to require a much lower amount 
of emission to originate from compact sources with a 
warmer dust emission peak. This result is not surpris- 
ing when one examines high angular resolution optical 
imaging of UGC 7321 and NGC 4244 (Figure IT]) in which 
dark clouds suggestive of molecular gas are observed. IC 
2233 shows some dark clouds in its central regions but 
less than the other two galaxies studied, reinforcing the 
suggestion that it is the galaxy poorest in interstellar 
molecular gas. 

From our total dust masses (Table [3| and the Hi 
and molecular hydrogen masses from theliterature (Ta- 
ble fTl) we can calculate global gas-to-dust ratios^ Af- 
ter accounting for the contribution of Helium (assuming 
Y = 0.25) we find global gas-to-dust ratios of ^ 666 
(UGC 7321), 734 (NGC 4244) and 1266 (IC 2233). These 
values are significantly higher than the Galactic value of 
~ 140 (|Whittet||1992|). The gas-to-dust ratios of UGC 



^'^ In the case of IC 2233 we use the upper limit H2 mass of 
1.4 X IO^Mq 



Large Magellanic Cloud ( Koornne ef 1982) and the Small 
Magellanic Cloud (SMC) ( [Bouc^t et al.|1985| while the 
value for IC 2233 is consistent with that ot the SMC. 

6.2. Model Limitations and Future Prospects 

It is somewhat unclear what effect the adoption of a 
clumpy dust and emissivity distribution w ould have on 
the FIR emi ssion of our galaxy sample. [Matthews fc| 
Wood ( 2001 ) found that in order to reproduce the opti- 
cal appearance of UGC 7321 it was necessary to allocate 
approximately 50% of the dust mass to a clumpy com- 
ponent. The effect that such a two phase diffuse dust 
distribution may have on the FIR emission would de- 
pend on the associated emissivity distributions. Such 
dense clumps could be "quiescent" and represent over- 
dense regions of the ISM illuminated only by the diffuse 
ISRF, but recent evidence suggests that such clouds are 
optically thin and will contribute little to th e attenua- 
tion and dust emission ( Molinari et al.pOlO l. However, 
if the dark clouds seen in H til ' imaging are star forming 
molecular clouds then they will be illuminated by both 
the ISRF and by young stars embedded within them. It 
is thought that in the case of LSB galaxies the optical 
depths of the clumps are similar to th ose found for diffuse 



molec ular clouds in the Milky Way ( Matthews & Wood 
2001|). In this case the clumpy component may provide 



a significant contribution to the FIR emission reducing 
the need to have a radially extended dust disk as is found 
using smooth density distribu tions. 



It has also been shown by Indebetouw et al. (20061 



that the observed SED of young stellar objects (YSUs 
that are embedded in clumpy circumstellar material can 
have a strong dependence on the observer's viewing an- 
gle (azimuth) and inclination. The observed MIR flux 
was found to vary by up to two orders of magnitude for 
the same highly embedded system viewed along a dif- 
ferent sightline. The effects are more severe at shorter 
wavelengths. Although the case of YSOs is very differ- 
ent to the galactic environment studied here it should be 
noted that the effects of a clumpy galactic ISM contain- 
ing optically thick clouds may cause viewing angle and 
inclination dependent effects on the MIR SEDs that are 
not accounted for in our smooth, axisymmetric models. 
Currently we have no observations of our galaxy sam- 
ple in the sub-mm wavelength regime. It has been found 
that in order to reproduce the sub-mm emission observed 
from some HSB galaxies i t becomes necessary to include 



et al. 2001 



an additional dust mass ( Popesc u et al.|2000 Misiriotis 
Popescu et aT 2011[ ). I'he additional dust 



mass reveals its presence through sub-mm emission that 
is underestimated by models which can adequately de- 
scribe the SED to ^ 100/im. Further evidence for an ex- 
tended cold dust component is also found by [Holwerda 



et al. ( 2005 ) based on attenuation of background galax- 



les m tace-on disks. Until observations become available 
we cannot be sure that our models do not lack cold dust 
emission at longer wavelengths. 

7. SUMMARY 

We have utilized multi-wavelength imaging and pho- 
tometry in conjunction with sophisticated Monte Carlo 
radiation transfer codes to investigate the structure of 
three edge-on, LSB disk galaxies. The galaxies have been 
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chosen to span a range in central, optical surface bright- 
ness and molecular hydrogen masses. 

We have been able to reproduce the global, optical 
appearance of all three galaxies using smooth emissivity 
and dust distributions. We find that the composition and 
size distribution of dust grains adopted, which are based 
on Milky Way extinction, provide a good match to the 
observed properties in our sample of LSB disk galaxies. 
Our models also reproduce the total emission at 70 and 
160/im for all three galaxies. However, the FIR morphol- 
ogy of our models appears more centrally concentrated 
than the more diffuse distribution suggested by the data. 
We find that the dust mass appears to be distributed in 
an exponential disk with a scale height comparable to 
or exceeding that of the stellar disk. This is in contrast 
to the findings for HSB galaxies where the dust disk is 
found have a ve rtical scale heig ht of approximately half 
the stellar disk (Xilouris et al. 1999). The comparable 
scale heights in tne dust and stellar disks is likely associ- 
ated with the increased stability of the ISM in LSB disks 



against vertical collapse ( Dalcanton et al.||2004 ) and the 
thin nature of the stellar disks, which suggests minimal 
dynamical heating. 

The dust masses and distributions derived suggest dust 
masses in the range 1.16 — 2.38 x 10^ Mq corresponding to 
face on, V band, optical depths between Tface = 0.034 — 
0.106. 

In future work we hope to develop our radiation trans- 
fer models to include small scale non-axisymmetric struc- 
tures which may shed further light on the structure of 
the ISM and star formation processes in LSB disk galax- 
ies. The inclusion of radial variations in the dust disk 
scale heights, associated with a flaring gas disk, could 
also prove important. Additional sub-mm observations 
may also allow us to uncover cold dust that is associated 
with the extended Hi, as has been found in som e HSB 
galaxies ( [Popescu fc Tuffs||2003[ |Hinz et al.||2006 1 . 
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